The dusty Nebula surrounding HR Car: a Spitzer view 
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ABSTRACT 

We present mid-IR observations of the Galactic Luminous Blue Variable 
(LBV) HR Car and its associated nebula carried out with the Spitzer Space 
Telescope using both IRAC and IRS, as part of a GTO program aimed to study 
stellar ejecta from evolved stars. Our observations reveal a rich mid-IR spectrum 
of the inner nebula showing both solid state and atomic gas signatures. Strong 
low-excitation atomic fine structure lines such as 26.0/im [Fe II] and 34.8/iin 
[Si II], indicate, for the first time, the presence of a PDR in this object class. 
While the physics and chemistry of the low-excitation gas appears to be domi- 
nated by photodissociation, a possible contribution due to shocks can be inferred 
from the evidence of gas phase Fe abundance enhancement. 

The presence of amorphous silicates, inferred from the observed characteristic 
broad feature at 10/xm located in the inner nebula, suggests that dust has formed 
during the LBV outburst. This is in contrast with the detection of crystalline dust 
in other probably more evolved Galactic LBVs, which is similar to the crystalline 
dust observed in red supergiants. This has been considered to be evidence of dust 
production during evolutionary phases prior to the outburst. 

Subject headings: circumstellar matter — infrared: stars — stars: early-type -stars: 
individual (HR CAR) — stars: winds, outflows 
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Introduction 



Luminous Blue Variables (LBVs) are luminous, massive stars, which represent a 
crucial and relatively short phase of massive star evolution between core-hydrogen burning 
O-type stars and helium burning Wolf-Rayet (W-R) stars. LBVs are quite rare objects 
i n our Galaxy - a recent census reported 12 confirmed members and 23 candidates 



( 1 Clark et al. 



20051 ) . From an observat ional point of view, LBVs are characterized by strong 



spectral and photometric variabili ty ( 



an increase of t 



re ma ss-loss rate (IDavidson 



( ILeitherer et al 



van Genderen 



200ll ). Such variability may indicate 



19871 ) or a change in the radius of the star 



19891 ) . Giant eruptions, during which a large amount of material is released 
(10~ 5 — lO _4 M ?/r _1 ), have been observed between prolonged periods of quiescence. Such 
events have been witnessed very rarely but the presence of extended circumstellar nebulae 
(LBVN) around Galactic and extra galactic LBVs suggests that they are a common aspect 



of LBV behavior (INota et al. 



19921 ) . Most of well known LBVNs contain large quantities of 



dust: they have been detected by IRAS, ISO, MSX a nd ground-based n e ar IR imaging as 



having excess IR and spatially resolved shells images (jTrams et al. 



1998 



Clark et al. 



20031 ). 



Egan et al 



2002 



According to recent evolutionary models, a massive star has to lose a huge quantity 
of its Main Sequence (MS) mass during its post-MS e yolution until i t reac hes its W-R 



Lamers et al. 



(1200 ll ) point out that 



phase during which it completely sheds its H envelope. 
LBVs and related transition objects (Ofpe/WN9 stars), which undergo massive mass-loss 
and eruptive events, may have a key role in the proposed scenario. There are, however, 
questions related to LBV phenomenon still under debate and, in particular, the nebula 
formation mechanism; i.e., the intensity, duration and geometry of single mass- loss event is 
not yet established. The full characterization of mass-loss properties during the LBV phase 
appears, therefore, crucial to assess their role in massive star evolution. 
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Much information on the LBV phenomenon can be derived by the study of the LBVNs, 
in particular the detailed analysis of gas and dust content in the LBVNs can be be used 
to tackle the problem of their origin. By studying the dust composition, it is possible to 
observe differences between LBVs that result from differences in the dust forming processes 
and/or the nebular abundances. The latter, as the dust chemistry would largely reflect the 
chemical composition of the gas phase, would help in localizing at which phase of stellar 
evolution the dust has been produced. Moreover, determining the distribution and the 
morphology of the gas and dust in the nebula may provide clues about the mass-loss history 
during and/or after the dust formation process. 

In order to begin to answer some of these questions, we have started a systematic study 
of a sample of LBVNs associated with Galactic LBVs using the Spitzer Space Telescope. 
The study is aimed to detect and resolve the faint dust shells that have been ejected from 
the central stars with mid-IR imaging, and to characterize the mineral composition of LBV 
ejecta via mid-IR spectra. In this paper we present results on the well-known Galactic LBV 
HR Car. 



1.1. HR Car 



HR Car is one of coolest Galac t ic LB V, classified as B2I. From a spectral analysis 
of the Balmer lines 



Manchado et al. 



(120021 ) have recently derived a stellar temperature of 
10,000-ft" and a luminosity of 5 x 1O 5 L , putting HR Car close to the 40M Q evolutionary 
track in the H-R diagram . The distance to HR Car (5 ± lkpc) has been derived by 



van Genderen et al 



(1l99ll ). using the reddening-distance method, which is consistent with 



th e value obtained from the k i nemat ic measurements of the Carina spiral arm (5.4±0.4/cpc) 



by 



Hutsemkers fc Van Prom 



(1991). 
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HR Car is surrounded by a faint, low-excitation nebula which is very difficult to observe 
in the optical due to the strong luminosity of the central object. Despite such difficulties, 
HR Ca r and its associate d neb ula have be e n ext ensively studied. Among the most recent 
results, 



Weis et al. 



(119971 ) and 



Nota et al. 



(119971 ). both based on optical high resolution 



coronographic imaging and spectroscopy, found that the nebula has a large-scale bipolar 
morphology, very reminiscent of the rj Carina nebula. Such lobes are expanding in the 



SE-NW direction (PA=135 degree) and have a diameter of ~ 18" ( INota et al. 



1993). On a 



smaller scale of a few arcsecs from the central object, a more compact, brighter nebula is 
evident, whose morp hological details, h owever, could not be derived from the coronographic 



optical (H a ) images. 



Nota et al. 



(119971 ) pointed out some degree of asymmetry, especially 
in the SE region of the inner nebula. The presence of some material close to the central 
object whose distribution i s not symmetric was also i nferred from the BVRI polarization 



measurements reported by iParthasarathy et al.l (120001 ) . 



Most of the information on the inner nebula can be obtained at longer wavelengths, 
where the difficulties due to the cont r ast be tween the bright central object and the faint 



nebula can be overcome. 



Voors et al. 



(119971 ) mapped the dust (broad-band 10 /zm) and the 
ionized gas (narrow-band [Ne II], 12.8 /im) of the nebula surrounding HR Car. The dusty 
nebula is quite compact, extending a few arcsecs from the central object, and reveals an 
asymmetric morphology. The ionized gas has the same extent as the dust, but a different 
distribution, implying a partial screening of the stellar ionization flux by the circumstellar 
dust. The most striking result is the difference between the morphology observed at smaller 
angular scale to that observed in the optical; this may imply that we are tracing different 
mass-loss episodes, which probab l y occ ur with a time-dependent geometry. The same 



conclusion was reached by IWhitd ()2000l ). who observed HR Car with the ATCA. The radio 
images reveal an extended ionized region whose morphology closely resembles that of the 
H a bipolar nebula plus an inner, compact nebula a few arcsecs across. The inner nebula is 
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strongly asymmetric and in complete disagreement with the large scale structure. Another 
puzzling feature of this LBV is that the central star doesn't appear to produce a sufficient 
quantity of ionizing photons to keep the nebula ionized, even assuming a previous period 
of higher temperature for the star. A possible explanation that takes in account both the 
observed degree of ioniza tion an d the s trong asymmetry of the inner nebula is a symbiotic 



model, as put forward by IWhitd (120001 ) . but the presence of a companion with the required 



properties has not yet been confirmed. 



Observations and Data Reduction 



The observ ations presented in this paper were carried out by usi ng the Infrared Ar ray 



Camera (IRAC; 



Fazio et al 



2004 ) a nd the Infrared Spe ctrograph (IRS; 



board the Spitzer Space Telescope (I Werner et al 



Houck et al. 



2004 ) on 



2004J). These observations were part of the 



GTO program (ID 30188) and were performed on 2006 July 10 (IRAC; AOR 17335040) and 
2006 June 28 (IRS; AOR 17335296). Both imaging and spectral observations were obtained 
at the central source position (RA 2 ooo = 10/i22m53.84s, DEC2000 = — 59d37m28.4s), which 
is determined, with good accuracy, from previous IR and radio imaging. 



2.1. IRAC 3.6, 4.5, 5.8 and 8.0 Images 



HR Car was observed with 8 IRAC frames per band using a medium cycling dither 
pattern at a frame time of 12 seconds, providing a field of view of approximately 6x6 arcmin 
in the final mosaics. The Basic Calib rated Data (BCD) p ipeline version S14.0.0 were 



used in the reduction, and IRACproc (jSchuster et al 



20061 ) was used to create the IRAC 



final images at a pixel scale of 0.6 arcsec/pixel (the original IRAC scale is ap proximately 



1.2 arcsec/pixel). The source structure revealed by previous mid-IR imaging (IVoors et al. 
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19971 ; iTrams et al.l 119971 ) is too compact to be well-resolved by IRAC. However, in these 
observations we are sensitive to fainter emission that is extended further from the star. 

At the 12sec frame time, the central star is saturated in all bands in the BCD. The 
PSF-subtraction technique in IRACproc was used to measure the flux in each of the bands. 
The PSF-subtracted mosaics are shown in Figure [TJ and the flux and magnitudes of the 
central star are given in Table [TJ The PSF-subtracted images have residuals, especially in 
the central regions, which are in part related to the slightly extended emission around the 
central star in the mid-IR. The photometry in Table Q] is the flux from the central star only 
and not the extended emission which is difficult to estimate close to the core. However, one 
can see evidence in the images for emission on the same spatial scale as seen in the Ha 



image (jNota et al 



19971 ). which is brighter at longer IRAC wavelengths. 



2.2. IRS Spectra 

The spectral observations were obtained by using the Short-High (SH) and Long-High 
(LH) modules (R ~ 600), which cover two spectral ranges going from 9.9 to 19.6 /im (SH) 
and from 18.7 to 37.2 /im (LH). For each module, three cycles were performed for each 
nod position in order to provide rejection of cosmic rays and other transient. From the 2D 
Basic Calibrated Data (IRS pipeline version 15.3.0) three spectra, for each nod position, 
were extracted and wavelength and flux calibrated using the Spitzer IRS Custom Extractor 
(SPICE) with the extended s ource extraction option. Using the Spitzer contributed software 



SMART flHigdon et al. 



2004 1. we cleaned the resulting ID spectra for residual bad pixels, 
spurious jumps, and glitches, and smoothed and merged them into one final spectrum per 
module. The resulting high resolution spectra for each module are shown in Fig. [5J 

The infrared emission from HR Car consists of narrow emission lines and a solid state 
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feature, at about 10 /im, superimposed on a thermal continuum. There is an expected 
mismatch between the SH and the LH spectra which is usually attributed to the difference 
in the two module apertures, as the SH has a slit size of 4.7 x 11.3" and the LH has a slit 
size of 11.1 x 22.3". The overplot of the IRS spectral slits for each high resolution module 
on the two Micron All Sky Survey (2MASS) K s (Fig. [3]) reveals that the central object is 
well contained in both apertures while only part of the inner nebula is contained in the 
SH slit. To quantify the correction for aperture losses, we can take advantage of the small 
spectral overlap between the two modules (from 18.7 to 19.6 /im) deriving a scaling factor 
of 1.87, at that wavelength range, to align the SH spectrum to the LH spectrum. 

However, we cannot a priori rule out a wavelength dependence of the missing flux, and 
a simple scale factor may be an oversimplification. We have therefore retrieved from the ISO 
Data Archive (IDA v9.0) a SWS spectrum (TDT=24900321), obtained over the full SWS 
grating range (2.38 — 45.2/xm), with a spectral resolution comparable to that of the IRS 
observations (~ 600), and with instrumental apertures ranging from 14" x 20" (2.4 — 12/mi) 
to 20" x 33" (29 — 45/im). Therefore, in the SWS observations the whole inner nebula w as 



Lamers et al. 



(1996) 



included in the aperture. This spectrum was already presented by 
but the data we retrieved from the archive are the result of further processing beyond the 
original pipeline with new, more refined algorithms (Highly Processed Data Products), 



providing a much improved spectrum with respect to the published one (ISloan et al. 



2003|). 



In Fig. H] the two Spitzer-IRS spectra are superimposed on the ISO data. A good 
agreement is evident up to ~ 10/xm. After that, there is a mismatch, with the amplitude 
increasing with wavelength, indicating that a simple scale factor cannot account for aperture 
losses. The IRS-LH is again in good agreement with the ISO spectrum up to ~ 27yum. At 
higher wavelengt hs in the ISO spectru m, there are additional apparent features which have 



been reported by 



Lamers et al. 



(119961 ). but these features are completely absent in the IRS 
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spectrum. We will further discuss the implications of these observed differences between 
the spectra in the following section. 



Solid state features 



The mid-IR spectrum of HR Car reveals a prominent feature that peaks at around 
10/im that is probably related to the 9.7/im feature due to the Si - str etching mode in 



amorphous silicates. This feature was also noted by 



Lamer s et al. 



(|l996l ). The presence of 



amorphous silicate in the stellar ejecta indicates an oxygen-rich (C/O < 1) environmen t, as 



we expect for LBV ejecta being enriched by CNO processed material (INota et al. 



19971 ). 



In the comparison of the IRS spectrum with the ISO SWS spectrum (FigHJ), we 
notice a reasonably good match between the two spectra at around lO/im. Since the two 
instruments have different apertures, the fact that in both spectra we see the same features 
with the same intensity implies that amorphous silicates must form very close to the central 
object, well inside the inner nebula, in a region whose extension must be comparable to or 
less than the IRS-SH aperture. At longer wavelenghts, we find that in the ISO spectrum 
there is a strong contribution beyond 27 /xm that is completely absent in the IRS spectrum. 
However the ISO spectrum is quite noisy in that wavelength range and it is difficult to 
assess the reliability of the observed broad structures. If we assume that the observed 
difference is related to the aperture, we may conclude that the emission features arise 
from a region located outside of the inner nebula, otherwise it should have been detected 
in the IRS spectrum. The features must be present in a region which is, at least partly, 
inside the are a seen by t he SWS aperture. T he nature of such a contribution is not clear. 



Lamers et al. 



(119961) and 



Waters et al. 



(119971 ) indicate a solid state origin, but no firm 



conclusion on the possible carrier/s was reached. 
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At this time, we can only speculate on the possible location of this material, but we 
cannot discriminate whether the spectral features are connected to some outer region of the 
HR Car nebula or to the ISM surrounding it. 



4. The emission line spectrum 



The mid-IR spectrum of HR Car shows many fine structure lines that, once the 
mechanisms responsible of the observed emission have been identified, may provide 
important constraints on the physical conditions in the line emitting region. The line-fitting 
routine in SMART was used to identify and measure the narrow lines. For each line a local 
continuum was defined by a single-order fit to the baseline and the line was then fitted 
using a single Gaussian. 

The most prominent lines are labeled in Fig[2j while the identification and the measured 
fluxes, including weaker lines (S/N < 20) are summarized in Table [2J Possible multiple 
identifications are also indicated. High ionization potential lines such as those of [S III] and 
[Fe III] can only originate in the H II region around the star, while the fine structure lines 
such as those of [Fe II] and [Si II] can form in a photodissociation region (PDR) surrounding 
the ionized part of the LBVN, where the gas is predominantly neutral but whose chemistry 



and thermal balance is contro 



13.6 eV (ITielens Hollenbach 



ed by FUV photons, i. e. with energies between 6 and 



19851 ). PDRs are bright in the far-IR dust continuum, in 



the PAH emission features, in some infrared fine-structure 



CO. CO emission detected in other LBV (INota et al. 



2002 



ines and in rotational lines of 



Rizzo et al. 



20081 ) supports the 



hypothesis of the presence of PDRs in the LBVNs. 
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4.1. [S III] lines 

The [S III] 33.48 frni and [S III] 18.71 fim lines can be used to derive the electron 
density in the ionized nebula associated with HR Car. Those lines originate from levels close 
in energy, and therefore their ratio will have almost no dependence on temperature, being 
only sensitive to their transition probabilities or collisional de-excitation, which depend on 
density. 

The [S III] 18.71 jum line is located in the overlap between the two high resolution 
modules. We can therefore safely apply the correction factor of 1.87 to take the aperture 
losses into account. From the observed ratio of the 33.48/im line flux to the 18.71/xm line 
flux, we derive, assuming a ele ctronic temperatur e of T e ~ 10 4 K, an electronic density of 



the H II region n e > 10 3 cm 3 (IHouck et al. 



1984J) . Because this line ratio is insensitive to 



temperature, the adopted value of T e will not affect the y a 



n e is in good agreement with those derived by 



Nota et al. 



ue of n P . The derived value of 



(119971 ). who obtained, from [S II] 
line ratios, a gradient in the electron density n e ranging from ~ 10 4 cm -3 in the center of 
the nebula to a value of ~ 10 3 cm -3 at a distance from the core of ~ 5". A similar value 
for the electron density (n e = 3200cm -3 ) was derived from the central core of the ionized 
nebula (6" x 8") by White (2000) from radio measurements. 

The agreement of our density determination with those derived by other authors in 
other spectral regions also supports the adopted scaling factor applied to the SH spectrum 
to account for aperture losses. The non-corrected 18.71/xm line flux would lead to a [S III] 
lines ratio that implies a very low density H II region (n e = 10 — 100cm -3 ), which is not 
consistent with results from visible and radio observations. 
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4.2. [Fe II] and [Si II] lines 

The [Fe II] 26.0/im and [Si II] 34.8/im lines may arise from the ionized region of the 
LBVN (H II region) and/or from the associated PDR. Therefore, we cannot a priori exclude 
a mixed contribution from these different regions to the observed lines intensities. The 
use of these lines as diagnostics of the physical condition i n the PDR is limited by the 



Kaufman et al. 



(J2006|) modeled 



knowledge of the real PDR contribution. In a recent paper, 
the intensities of typical infrared fine-structure lines that may originate from both the H II 
region and the associated PDR in order to evaluate the contribution to such diagnostics 
from the two different regions. They concluded that, in the hypothesis of thermal pressure 
balance between the H II regions and the PDR and normal metallicity, the fine structure 
line of [C II] at 158.0/im is always dominated by the PDR, while the [Fe II] 26.0/im and 
[Si II] 34.8/im emission is dominated by the PDR if the H II region has a h igh electron 



density n e > 10cm 3 . A similar conclusion was reached by 



Abel et al. 



fl2005h . 



Assuming the n e for the H II region associated with HR Car as derived from the 
[S III] line ratio, we conclude that the [Fe II] 26.0/im and [Si II] 34.8/im lines observed in 
the mid-IR spectrum of HR Car are PDR-dominated and can be safely used to derive the 
physical conditions of this region. For this purpose we need to compare the observed line 
intensities with those predicted by the theoretical PDR models. 

The physics and chemistry of a PDR are governed by the incident FUV (6eV < hv < 
13.6eV) radiation, often expressed as relative to the background interstellar radiation field 
assumed to be 1.6 x 10~ 3 erg cm" 2 s _1 : 

° I\d\ 

q _ J 912A A 

1.6 x 10~ 3 (erg cm~ 2 s _1 ) 

and by the total density n a (cm -3 ). 

In Figures [5] and [6], the line intensities for [Si II] 34.8/im and [Fe II] 26.0/im computed 



13 



with the model of Kaufman et al.l (120061 ) are shown as function of G a for several values of 
n Q . The measured fluxes of the [Fe II] 26.0/im and [Si II] 34.8/mi lines correspond to a line 
intensity of 1.04 x 10~ 3 (erg cm~ 2 s~ x sr~ l ) and 1.42 x 10 _3 (erg cmT 2 s^ 1 sr^ 1 ) respectively, 
if we assume that the PDR completely fills the LH aperture (Qpdr — 0.57 x 10~ 8 sr). The 
shaded regions in each figure indicate our line measurement with the associated uncertainty. 

Different combinations of the space parameters (G Q , n a ) can explain the observed 
line intensity of [Si II] at 34.8yum, going from a high density (n a ~ 10 6 ) and a moderate 
incident FUV radiation log(G Q ) ~ 3.5 to a low density (n ~ 10 3 ' 5 ) and a high incident 
FUV radiation log(G ) ~ 6 (Figure E]). On the other hand, to explain the observed line 
intensity of [Fe II] at 26.0/im, a very high density (n > 10 7 ) and a high incident FUV 
radiation (log(G ) > 4.5) are necessary. 



Because the [Fe I 



(ITielens fe Hollenbach 



1 line at 26.0/xm and the [Si II] line at 34.8/zm are optically thin 



19851 ). the modeled line intensities scale linearly with the adopted 



gas phase abundance. This means that to have physical conditions of the PDR consistent 
with both the observed line intensities, the gas phase F e abun dance in the nebula should be 



Kaufman et al. 



(120061 ) standard model typical for 



higher than the values adopted by the 
the diffuse ISM (1.7 x lO" 7 ). 

It is quite difficult to quantify the gas phase Fe abundance necessary to reproduce 
th e observed 26.0/xm [Fe II] lin e intensity without detailed modelling. The PDR model 



of 



Tielens fe Hollenbach! (119851 ). which assumes a gas phase Fe abundance of 2.5 x 10 7 , 



foresees a 26.0/xm [Fe II] line intensity of 4 x 10~ 4 erg cm _2 s _1 sr~ 1 , for G a = 10 5 and 
n = 10 4 " 5 . Therefore, with an increment of ~ 50% of the gas phase Fe abundance, the 
expected 26.0/im [Fe II] line intensity is enhanced of a factor of four, shifting the model 
curves shown in Fig [6] towards our observed values. 
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5. Conclusions 

We have presented mid-IR observations of the Galactic LBV HR Car and its associated 
nebula carried out with Spitzer using both IRAC and IRS, as part of a GTO program 
aimed to study stellar ejecta from evolved stars. 

The central object in HR Car is too bright in the mid-IR and the warmer inner 
dusty nebula too compact to get reliable mapping with IRAC; still, the adopted observing 
strategy has allowed us to observe the presence of faint structures, brighter at longer IRAC 
wavelengths, whose extension is comparable to those of inner ionized nebula traced by 
H a and the radio. The more extended nebula, which, in the mid-IR, can be detected 
only through its free-free continuum (no strong recombination lines are evident in the ISO 
spectrum for A < 10/im), is expected to be too faint to be observable at the integration 
time used. 

We have further obtained a mid-IR spectrum of the inner nebula and a few arcsecs 
across the central object. The spectrum reveals quite a number of emission lines, some of 
which are formed in the ionized part of the nebula. Particularly interesting is the presence 
of low-excitation atomic fine structure lines such as 26.0yum [Fe II] and 34.8//m [Si II], lines 
whose detection may indicate the presence of a PDR. The comparison between the observed 
line intensities and those derived from recent PDR models supports the existence of a PDR 
surrounding the ionized fraction of the nebula. However, to obtain a PDR with physical 
characteristics that are able to reproduce both [Fe II] and [Si II] line intensities, a Fe gas 
phase abundance higher than that of the ISM is necessary. 

In the ISM, Fe is usually highly depleted from the gas phase because of condensation 
onto dust grains. In this context, a higher Fe gas phase abundance with respect to that of 
the ISM may be considered as indirect evidence of shocks occurring in the nebula resulting 
in grain processing that release Fe atoms into the gas phase. 
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Such an enhanced Fe abundance has been observed in the radio Arc Bubble of the 
Galactic Center and has been interpreted in terms of additional Fe release d to the gas 



phase as a consequence of grain destruction by shocks (jSimpson et al. 



20071 ). Such shocks 



are probably produced in the strong stellar winds of the Quintuplet Cluste r stars , which 



Smith! (120021 ) 



also includes a LBV, the Pistol star. Shock-heating has been proposed by 
to explain the bright 1.643/iin [Fe II] line observed in a small sample of LBVs. While 
the presence of such a line seems to be a common property of LBV nebulae, its possible 
excitation mechanism has not been clearly established. We conclude that in the case 
of HR Car, the photodissociation dominates the physics and the chemistry of its outer 
envelope, but shocks also play an important role as inferred from the enhanced gas phase 
Fe abundance. 

The analysis of the IRS spectrum has allowed us to derive some clues on the mineralogy 
and spatial distribution of the dust. There is a strong signature of amorphous silicate, 
which appears to be localized within ~ 10" (0.3pc, assuming a distance of 5.4kpc) from 
the central object, in a region spatially coincident with the optical, IR and radio inner 
nebula. Our result is consistent with the amorph o us sili cates being localized in the compact 



3") nebula mapped at 10 jum by 



Voors et al. 



(119971 ). Such a strong emission band is 



probably also responsible for the diffuse emission detected in the 8.0 ± 1.5/xm IRAC band. 



The spatia 



Nota et al 



coinc idenc e among th e IRAC structures with the inner H a nebula reported by 



(I19971 ) and 



Wcis et al 



(119971 ) strongly supports our conclusion that ionized gas 



and amorphous silicate are localized within the same area. 

Relative to the region of the nebula covered by the IRS apertures, we cannot confirm 
the presence of crystalline features indicated by previous ISO observations. The presence of 
amorphous silicates, together with the lack of crystalline silicates in the inner part of the 
nebula, suggest that there has been recent dust formation in HR Car that has condensed 
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during the LBV eruptions. This is in contrast with other Galactic LBVs, such as AR Car 
and WR751, where crystalline dust has been observed. The similarity of the crystalline 
dust in those objects to the dust observed in red supergiants has been considered evidence 
of dust pro duction in a previo us evolutionary phase of these massive stars, prior to the LBV 



outbursts (jWaters et al. 



1998|) 
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Fig. 1. — The PSF-subtracted images of HR Car based on the 12 sec frames obtained in 
each of the four IRAC channels. The core of the source was saturated and therefore appears 
as a black region in the center of the image. The images were rotated approximately 123.3 
degrees counterclockwise to align the m to R.A. and De c, shown, therefore the array banding, 



muxbleed, and bandwidth artifacts (jHora et al. 



20041 ) are also rotated from their nominal 



BCD orientation and appea r as dark or 



frame. The Ha image from 



Nota et al 



brigh t residuals that are roughly diagonal in the 



(jl997l ) is shown to the right of the IRAC images. 



There is low-level extended emission is visible in the images which is similar to that seen in 
the Ha image. The residual rings near the core result from a mismatch between the PSF 
and HR Car images, possibly indicating that the core is slightly broader than a point source. 
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Fig. 2. — IRS High-resolution spectrum of HR Car. See the electronic edition of the Journal 
for a color version of this figure. 
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Fig. 3. — Overplot of the IRS spectral slits, for both high resolution modules, on the two 
Micron All Sky Survey (2MASS) K s image. The smaller aperture (blue box in the electronic 
edition) corresponds to the SH aperture, the larger one (red box in the electronic edition) is 
the LH aperture. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 4. — The IRS/SH and LH spectra superimposed to the ISO-SWS01 archive spectrum. 
See the electronic edition of the Journal for a color version of this figure. 
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Fig. 5. — The line intensity for [Si II] 34.8/xm computed with the model of iKaufman et al. 



20061 1 is shown as function of G Q for several values of n . The shaded region indicates our 



line measurement with associated uncertainty. 
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Fig. 6. — Tha same as fig. [5] but for [Fe II] 26.0/iin line. 
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Table 1: HR Car Fluxes measured in the IRAC channels 



IRAC Band Flux Density (Jy) Mag 

3.6 0.985 6.14 

4.5 0.848 5.82 

5.8 0.714 5.51 

8.0 0.910 4.62 
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Table 2: Line Identifications and Measured Fluxes 
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a Obscrvcd fluxes in units of 10~ 14 erg cm~ 2 s 

h Flux uncertainties < 10% 

Tlux uncertainties between 10% and 20% 



